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ABSTRACT 

We report the detection of linear polarization from Sgr A* at 750, 850, 1350 and 
2000/um which confirms the contribution of synchrotron radiation. From the lack of po- 
larization at longer wavelengths it appears to arise in the millimetre/sub-millimetre ex- 
cess. There are large position angle changes between the millimetre and sub-millimetre 
results and these are discussed in terms of a polarized dust contribution in the sub- 
millimetre and various synchrotron models. In the model which best explains the data 
the synchrotron radiation from the excess is self-absorbed in the millimetre region and 
becomes optically thin in the sub-millimetre. This implies that the excess arises in an 
extremely compact source ~ 2 Schwarzschild radii. 

Subject headings: Galaxy: centre-Synchrotron radiation-Sgr A* 

1. Introduction 

Observations of stellar proper motions in the vicinity of Sgr A* (Ghez et al. 1998), the non- 
thermal radio source at the apparent centre of the Galaxy, show that its mass of ~ 2.5 x 10 6 M@ 
is highly compact on a scale < O.Olpc, reinforcing its claim to be a massive black-hole candidate. 
The spectrum of Sgr A* extends as a rough power law ~ v a , where a lies between j and |, from a 
low-frequency turn-over at a few GHz up to ~ 100GHz, (e.g. Mezger, Duschl and Zylka 1996 and 
references therein), and above this frequency there is evidence of a mm to sub-mm excess over this 
power law extending almost to the atmospheric cut-off near 1000 GHz (Serabyn et al. 1997, Falcke 
et al. 1998). The nature of the mm/sub-mm excess has been discussed by Serabyn et al. (1997) 
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and the possibility that it is due to dust is effectively eliminated: at 1.3mm its size is less than 1 
arcsecond with an implied brightness temperature in excess of 100 K, inconsistent with the shorter 
wavelength data, and the implication is that there are two separate synchrotron components. 

The radio spectrum of Sgr A* has for some years been modelled in various ways in terms 
of synchrotron radiation but searches for linear polarization, the characteristic signature of this 
mechanism, have only recently been reported. Bower et al. (1999a,b) used the VLA at 4.8, 8.4, 
22 and 43 GHz in its spectropolarimetric mode and the BIMA array at 86 GHz finding upper 
limits of <0.2% at the lower frequencies and <1% at 86 GHz. These small upper limits on linear 
polarization, given the sensitivity of the observations to large rotation measures due to Faraday 
rotation (RM ~ 10 7 rad m -2 at 8.4 GHz), are difficult to account for and this is discussed at length 
by Bower et al. (1999a). Apart from this there has been only a marginal detection at 800 fim of 
4.9 ± 3.2% at 129 ±19° (Flett and Murray 1991) and a report of circular polarization at 4.8 GHz 
by Bower et al. (1999c). 

In mm to sub-mm polarimetric imaging studies of the central 15 pc of the Galaxy we have 
detected linear polarization in Sgr A* at 750, 850, 1350 and 2000 /im which we present here. 
Discussion of the field distribution in the circum-nuclear disk (CND) and neighbouring molecular 
clouds as revealed by these observations will be presented separately. 

2. Observations and Results 

The observations were made using the SCUBA camera (Holland et al. 1999) and its polarimeter 
(Greaves et al. 1999) on the 15m James Clerk Maxwell Telescope (JCMT) in Hawaii. SCUBA 
observes simultaneously with arrays of 37 and 91 bolometers, at either 850 and 450 /im respectively, 
or 750 and 350 /im with a change of filters. Polarimetric observations were made in imaging mode 
at all four submillimetre wavelengths, and single-point polarimetry towards Sgr A* was also done 
at 1.35 and 2mm, using single bolometers and supplementing the data with small grid maps. The 
arrays have fields of view of 2.3 arcmin, and full-width half-maximum beam sizes range from 7-8" 
at the shortest wavelengths to 34" at 2 mm (Table 1). 

Observations were made by chopping the secondary mirror at 7.8 Hz to remove the mean sky 
level, and nodding between left and right beams at slower rates to take out sky gradients. For 
polarimetry, a rotating half-wave plate and fixed etched grid were used to modulate the signal 
seen by the detectors. Polarimetry at 750 and 850 /an resulted in a Nyquist-sampled image after 
32s of integration, after which the waveplate was stepped by 22.5°. The single-pixel data were 
observed in a similar manner but with 8s of integration per waveplate angle. Total integration 
times were 5-10 complete waveplate cycles, or 15-85 minutes. The data were reduced using the 
SURF (Jenness and Lightfoot 1998) and POLPACK (Starlink user note 223) software packages; 
for each pixel a sinusoidal modulation is fitted to deduce the source percentage and direction (p, 6) 
of the polarization. 
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From observations of Saturn and Uranus, during the run and from archival data, the instru- 
mental polarization is found to be close to 1% and known to ±0.1-0.2 % at a given wavelength and 
elevation. At 450 and 1350/mi the instrumental polarizations are larger (3.5 and 2.7% respectively 
due to wind blind structure) but still known to 0.2 to 0.8% respectively. The planets were assumed 
to be unpolarized and a limited check was made on this for Saturn where sky rotation can be used 
to separate instrumental from any planetary polarization. While we cannot exclude the possibility 
that Saturn is slightly polarized this introduces an error of less than 0.3% in the instrumental 
polarization subtraction. The instrumental polarization changes by only 0.1% over the elevation 
range used and the uncertainties in the corrected polarization is considered good to ±0.1 — 0.2%. 
Errors in the absolute orientation of the waveplate are one degree or less. 

Intensity calibration at 1350 and 2000/mi was obtained using Uranus taking the brightness 
temperatures as 96K and 110K respectively, and at 850 fim from Mars taking Tg = 208K. Since 
the emphasis is on the polarimetric results, we did not perform extremely accurate calibration 
observations, but typical uncertainties are only about ±10 %. For the 750 fim data Saturn was 
used for calibration taking Tg = 123K and a greater error was introduced since this planet is larger 
than the beam; we note that the 750 /an flux for Sgr A* appears anomalously low. Atmospheric 
tranmission was measured with skydips. 

A fundamental limitation is the maximum chop throw of 180" (in this case at a position angle 
of 145° east of north), which proved to be insufficient at 350 and 450 fim. Surrounding dust 
emission has a steep spectral index (Pierce-Price et al.. 2000) resulting in increasing off-beam 
contamination at these shortest wavelengths where Sgr A* was not detected above background. At 
850 ^m the off-beam signals are only 15% of the Sgr A* flux and these uncertainties should be 
similar or smaller at 750, 1350 and 2000 /Ltm at all of which Sgr A* was clearly detected. During 
the observation period the field of view rotates but the chop orientation on the sky is maintained; 
thus every spatial point in the map retains the same chop positions in a and 5 although different 
array bolometers may be involved. The results of polarimetric imaging in the vicinity of Sgr A* 
are presented in Table 1 and the total flux at the position of Sgr A* is given in column 3. 

There will in general be two contributions to the background, namely from cool dust in the 
inner cavity and the edges of the CND and free-free emission from the HII region Sgr A West, 
predominantly from the central ionized filaments. 

Dust and free-free emission contribute significantly to the observed flux but in this wavelength 
region there is no independent evidence at sufficient resolution to determine the contribution of 
the dust emission component in the central beams. In the sub-mm the flux in annuli from one 
to two beam radii were determined from the images and at 1350 and 2000/um observations in 
eight beams circumferentially distributed about Sgr A* sampled the ambient flux between 1.7— 
3.7 beam radii from Sgr A*. Estimates of the free- free contribution in the central beams and 
these peripheral areas were made from the 2 arcsecond resolution 3.6cm radio continuum maps of 
Roberts and Goss (1993), after removing the Sgr A* point source itself. These fluxes were then 



-4 - 



scaled to the present wavelengths as v , appropriate for optically thin free-free emission, and the 
approximation is made that the central dust contribution is just that derived from the peripheral 
regions after subtraction of the free- free component (column 7 in the table). The flux from SgrA* 
after subtraction of the dust and free-free contributions (column 8) will also include any local dust 
excess (or deficit) over background. Since at 2mm the correction for free-free is large due to the 
large beam size the reliability of the non-thermal flux estimate at this wavelength will be affected. 

The observed E vector polarizations and position angles (north through east) are listed in 
Table 1 in columns 9 and 10. At 450/mi there is general polarization over the central 30 arcseconds 
with no significant change at the position of Sgr A* and its average ~ 3% at a position angle of 
100° (E vector) is attributed to dust. At wavelengths between 750 and 2000/mi the polarization 
at the position of Sgr A* is well detected at the lO-u level. In the sub-mm it stands out from its 
surroundings and is clearly associated with the flux peak at the position of Sgr A*. This can be 
seen in Fig 1 which shows the central 80 x 100" region about Sgr A* at 850/Lim. At 1350 and 
2000/mi the observations are from single pointed observations where the HII region Sgr A West 
and Sgr A* are the dominant sources in the beam and Sgr A* is the only significant contributor to 
polarization. 

3. Interpretation 

At 2000/mi, and to a lesser extent at 1350/xm, there are only two contributors to the flux: 
diffuse free-free emission and Sgr A* itself. Since the former is unpolarized the polarization observed 
confirms the synchrotron nature of Sgr A*; it also suggests that the polarization is only associated 
with the mm/sub-mm excess since it has not been observed at wavelengths of 3.5mm and longer 
where the power-law spectrum dominates. Note that at 2000^m the largest contributer is from 
free-free emission and therefore the intrinsic polarization of SgrA* must be large. 

A large difference of position angle between the mm and sub-mm observations is evident in 
Table 1 and various ways of understanding this are next investigated. 

Although the mm position angles themselves are very close the three shortest wavelength 
position angles define a straight line when plotted against A 2 ; if a rotation of tt is subtracted from 
the 2mm angle a rough linear relation is maintained. However the rotation measure then indicated 
is so large (—1.44 x 10 6 rad/m 2 = vr/(A 2 — A 2 35 )) that polarization within the band pass of 40GHz 
at 2mm would be undetectable. Therefore the closeness of the mm position angles indeed indicate 
that there is little if any Faraday rotation. 

Changes in polarization and position angle with time are well documented in the compact 
cores of blazars and some AGN (e.g. Saikia and Salter 1988) and also have been observed at mm 
wavelengths (Nartallo et al. 1998). Sgr A* is known to be variable on a time scale of months. 
While there is a separation of some months between the observations presented here, those at 750, 
1350 and 2000/im were all made in August 1999 within a week and these do show the large position 
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angle difference between the mm and sub-mm polarizations. Furthermore, the 850//m observation 
in March also lies on the same trend of position angle with wavelength. Variability is therefore 
unlikely to have been a significant influence on the position angle, although it seems possible that 
Sgr A* may have increased in brightness during this period. 

Since the sub-mm dust emission is polarized it will affect the observed polarization and it is 
tempting to try to attribute the whole position angle change to dust. 

Denoting by I Q , Id and I s the observed central flux, the estimated dust and Sgr A* components 
respectively and by q a , qd and q s their Stokes q parameters, the synchrotron component q s can be 
found from 

IoQo = hQd + I s q s 
and similarly for the Stokes u components. 

It is difficult to estimate the dust contribution to polarization from the 750 and 850^m obser- 
vations because of the proximity of Sgr A* and its high intrinsic polarization. Instead we use the 
450/Ltm observations in which Sgr A* is not detected, and the polarization is ~ 3% at a position 
angle of ~ 100°. There are some problems here with polarized flux in the reference beams but other 
indicators of dust alignment in the central parsec are consistent with this estimate: observations at 
350/um (Novak et al. 2000) yield 1.9% at 87° near the position of SgrA*, and at lOO^rn the results 
are similar (Hildebrand et al. 1993) though both these are in a larger beam; the only information 
at high resolution is from the mid-infrared (Aitken et al. 1996) where integrations over a central 
14 arcsecond diameter yield 2.4% at 125°, admittedly weighted to warm dust. 

The intrinsic polarization resulting from removing the unpolarized free-free contribution and 
adopting 3% at 100° for dust emission is shown in Table 1 columns 11 and 12. At mm wavelengths 
these numbers are insensitive to the uncertainties in the adopted dust polarization and even in 
the submm the range of dust polarizations given in the previous paragraph changes the intrinsic 
polarization by ±3% at most and less than ±3°. Similarly uncertainties in the relative fractions of 
the three flux components has little effect on the intrinsic position angle but changing dilution does 
introduce errors to the intrinsic polarization fraction and at 750 and 2000//m the upper bounds are 
poorly defined. In Table 1 the errors to the intrinsic polarization are derived by assuming ±15% 
errors to the free-free and ±20% errors to the dust contributions except at 2000/mi where the 
adopted error is ±40% (these dust contribution errors are derived from the observed variation in 
the peripheral flux). It is clear that dust polarization cannot explain the abrupt change of position 
angle between the mm and sub-mm, rather it is increased and it approaches ir/2. 

We can turn the above argument around to determine what dust polarization is needed in the 
sub-mm to cause the position angle shift given the intrinsic synchrotron polarization of the mm 
results, say ~ 10% at 85°. With the 850^m fluxes and polarization given in Table 1 the required 
dust polarization is then 12.0% at 163°, larger than any dust emission polarization yet observed. 
It might be argued that this could be due to unusual conditions close to Sgr A* and an additional 
dust component in the central beam. Any such emission is already included in the Table 1 entry 
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for Sgr A* and because the same polarized intensity (~ 450m Jy) has to be supplied from much less 
dust emission the required polarization would be in excess of 25%. 

An assumption made so far is that throughout the wavelength region of these observations the 
synchrotron radiation remains either optically thin or self-absorbed. Most models of the synchrotron 
SED in Sgr A* consider the mm/sub-mm excess to be self absorbed but it may well be that at the 
shortest wavelengths it becomes optically thin. In that case the synchrotron position angle shifts 
through 90° and becomes close to 175° in the sub-mm. The dust corrected sub-mm points are now 
discrepant by 6 to 14° and this could be attributed to the crudeness of the correction. 

Approximating the synchrotron emitting region by a simple slab with constant energy density 
and magnetic field we find (e.g. Ginzburg and Syrovatskii 1965, 1969, Pacholczyk and Swihart 
1967) that for the polarization to change sign near 1mm then u ssa , the synchrotron self absorption 
frequency, lies in the range 0.6-0. 7mm, close to the peak of the mm/sub-mm excess, and the 
required fields, densities and energies are similar to those in the compact source models of Beckert 
and Duschl (1997). The polarization transition of this simple model is quite sharp, as is observed, 
with the polarization quickly approaching saturation on the long wavelength side, changing sign 
near 1mm and steadily rising at shorter wavelengths. 

A variant on this explanation is that self-absorption results in different optical depths with 
possibly different field distributions being probed as a function of wavelength. In that case, however, 
one should encounter diminishing polarization towards shorter wavelengths due to the superposition 
of a range of field orientations, and this is not observed. 

4. Discussion 

Synchrotron radiation from Sgr A* has been detected between 750//m and 2000/um with an 
intrinsic polarization ~ 10%. This polarization arises from the spectral region of the mm/sub-mm 
excess. The present observations are compounded by the presence of polarized emission from dust 
and dilution by free-free emission in these large beams. Although this has led to a number of 
possible explanations for the large position angle shift between the mm and sub-mm results all but 
one of these can be eliminated: the closeness of the position angles at 1.35 and 2mm require that 
there is little or no Faraday rotation, dust polarization can only produce a small shift which is in 
the wrong sense, and the shift is not due to variability. Only a transition between optically thin 
synchrotron radiation in the sub-mm and self-absorption in the mm region appears as a plausible 
explanation of the shift. Future small beam observations will be needed to confirm this result. 

A puzzle re-emphasized here is the lack of polarization at 3.5mm (Bower et al.. 1999a, b) and 
longer wavelengths in view of the large intrinsic polarization of Sgr A* at 1.35 and 2mm. Sufficient 
dilution of the mm/sub-mm excess to bring it below 1% requires not only a steeply falling spectrum 
with decreasing frequency for the excess, where i/ 5 / 2 is expected for self absorption, but also that 
the cut-off frequency of the power-law section is below 2mm to ensure a fast rising flux of the 
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power law region. It is possible, as well, that a flare in Sgr A* at 100GHz in March 1998 (Tsuboi, 
Miyazaki and Tsutsumi 1998) has affected the result of Bower et al. 1999b. 

The synchrotron self-absorption transition in the sub-mm proposed here is consistent with the 
model of the excess by Beckert and Duschl (1997), and requires the excess to arise from a compact 
source ~ 2 Schwarzschild radii in size. Also because of the transition the ir/2 ambiguity of the 
relationship between position angle and magnetic field is removed and hence the field orientation 
is close to east-west, since at mm wavelengths the emission is considered self-absorbed. 

5. Conclusions 

We report the observation of millimetre and sub-millimetre polarization from Sgr A*, confirm- 
ing the role of synchrotron radiation. The polarization is a property of the mm/sub-mm excess, 
demonstrating that the excess is real and not an artefact of variability and that it and the power- 
law spectrum arise in distinct structures. There is a large position angle shift between the mm 
and sub-mm observations. This can be explained by a transition between optically thin and self 
absorbed synchrotron radiation near 1 mm. Such a high self-absorption frequency implies a very 
compact source ~ 2Rs- 
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Fig. 1. — Greyscale image of 3.6 cm emission (Roberts & Goss 1993, courtesy of ADIL) overlaid 
with contours of the 850 /xm emission from the central region around Sgr A*, which is clearly seen 
as a point source in the centre of this image. The short lines on the figure show the magnitude and 
direction of the sub-millimetre polarisation E-vector. Note how the magnitude of the polarisation 
is larger at the position of Sgr A* (more than one large vector is seen as the beam is bigger than the 
vector spacing in this figure). Coordinates for Sgr A* are 17 h 45 m 40.049 s , -29° 00' 27.98" (J2000). 
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Table 1. Sgr A* fluxes 



A beam central flux (Jy/beam) ambient central Sgr A* obs. polarsn. intr. polarsn. date 

fim arcsec total excess free-free free-free dust P(%) P(%) $(°) 

1 2 3 4 a 5 6 7 b 8 9 10 11 12 13 

=3-4-6 =4+6-5 



750 


12. 5 C 


~5 


~2.5 


2.5 


1 


~1.5 


~ 1 ±0.5 


3.8±0.4 


164±3 


22±f 


169±3 


27 Aug 99 


850 


14 


8.2 


3.6 


2.25 


0.95 


3.65 


2.3±0.9 


3.3±0.2 


150±2 




161±3 


25 Mar 99 


1350 


22.5 


6.5 


4.9 


3.55 


0.85 


0.75 


2.2±0.5 


4.1±0.4 


89±3 


nil 


88±3 


22 Aug 99 


2000 


33.5 


8.0 


6.4 


5.8 


1.2 


0.4 


1.8±0.9 


2.9±0.3 


84±3 


12±» 


83±3 


22 Aug 99 



a (ccntral excess) = (central total)- (ambient) 

b flux from the central dust component is taken as equal to the flux from the ambient dust 
c fluxes were measured in a 16 arcsec beam 

Note. — Errors on observed P and 6 are la statistical errors. See text for uncertainty estimates on Sgr A* flux and intrinsic 
polarisations. 



